Typeset with jpsj2.cls <ver.l.2.2b> 



Letter 



00 
O 

o 

(N 

> 
O 

:^ 

in 

(N 



Magnetism and Superconductivity in the Two-Dimensional 16 Band d-p Model for 

Iron-Based Superconductors 

Yuki Yanagi*, Youichi Yamakawa and Yoshiaki Ono 

Department of Physics, Niigata University, Ikarashi, Niigata 950-2181, Japan 

The electronic states of the Fe2As2 plane in iron-based superconductors are investigated on 
the basis of the two-dimensional 16-band d-p model which includes the Coulomb interaction 
on a Fe site: the intra- and inter-orbital direct terms U and U' , the Hund's coupling J and 
the pair- transfer J'. Using the random phase approximation (RPA), we obtain the magnetic 
phase diagram including the stripe and the incommensurate order on the U'-J plane. We also 
solve the superconducting gap equation within the RPA and find that, for large J, the most 
favorable pairing symmetry is extended s-wave whose order parameter changes its sign between 
the hole pockets and the electron pockets, while it is dxy-wme for small J. 
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The newly discovered iron-based superconduc- 
tors^- 2) RFePnOi-xF^ (R=Rare Earth, Pn=As, P) with 
a transition temperature up to Tc = 55K^^ have at- 
tracted much attention. The F nondoped samples exhibit 
the stripe-type antiferromagnetic order with a transition 
temperature 134K and a magnetic moment ~ 
at low temperature. With increasing F doping, the sys- 
tem becomes metallic and the antiferromagnetic order 
disappears, and then, the superconductivity emerges 
for X ~ 0.11 with Tc ^ 26K. Specific features of 
the systems are two-dimensionality of the conducting 
Fe2As2 plane and the orbital degrees of freedom in Fe^+ 
(3d^).^'^^ The pairing symmetry together with the mech- 
anism of the superconductivity is one of the most signif- 
icant issues. 

The NMR Knight shift measurements revealed that 
the superconductivity of the systems is the spin-singlet 
pairing.^' ^) Fully gapped superconducting states have 
been predicted by various experiments such as the pen- 
etration depth, '^^ the specific heat,®) the angle resolved 
photoemission spectroscopy (ARPES)^"^^) and the im- 
purity effect on Tc.^' In contrast to the above mentioned 
experiments, the NMR relaxation rate shows the power 
low behavior l/Ti oc below Tc,^^) suggesting the su- 
perconducting gaps with line nodes. 

Theoretically, the first principle calculations have 
predicted that the nondoped system is metallic with 
two or three concentric hole Fermi surfaces around P 
(k = (0,0)) point and two elliptical electron Fermi sur- 
faces around M (k = (tTjTt)) point. Mazin et al. 
suggested that the spin-singlet extended s-wave pairing 
whose order parameter changes its sign between the hole 
pockets and the electron pockets is favored due to the an- 
tiferromagnetic spin fluctuations."'^^' ^^-^ According to the 
weak coupling approaches based on multi-orbital Hub- 



bard models, the extended s-wave pairing or the 
dxy-vfa.ve pairing is expected to emerge. The details of 
the band structure and the Fermi surface are crucial for 
determining the pairing symmetry. Therefore, theoretical 
studies based on a more realistic model which includes 
both the Fe 3d orbitals and the As 4p orbitals, so called 
d-p model, are highly desired. 

In the previous paper, ^"^^ we have investigated the 
pairing symmetry of the two-dimensional 16-band d-p 
model by using the random phase approximation (RPA) . 
It has been found that, for a larger value of J/U', the 
most favorable paring symmetry is extended s-wave, 
while, for a smaller value of J/U' , it is dxy-'wave. How- 
ever, the detailed electronic states in the whole param- 
eter region of U' and J have not been discussed there. 
The purpose of this paper is to obtain the detailed phase 
diagram including the magnetism and the superconduc- 
tivity in the U'-J plane. 

First of all, we perform the density functional cal- 
culation for LaFeAsO with the generalized gradient ap- 
proximation of Perdew, Burke and Ernzerhof by using 
the WIEN2k package, where the lattice parameters 
(a = 4.03268A, c = 8.7411lA) and the internal coordi- 
nates {zLa — 0.14134, ZAs = 0.65166) are experimentally 
determined. Considering that there are two distinct Fe 
and As sites in the crystallographic unit cell, we then de- 
rive the two-dimensional 16-band d-p model, where 
3d orbitals {d3z^_r2, dx2-y2, dxy, dyz, dzx) of two Fe 
atoms (Fe^=A, Ve^=B) and 4p orbitals {px, Py, Pz) of 
two As atoms are explicitly included. We note that x,y 
axes are rotated by 45 degrees from the direction along 
Fe-Fe bonds. The noninteracting part of the d-p model 
is given by the following tight-binding Hamiltonian, 



H. = 



tlrjdila- 
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where d^a is the annihilation operator for Fe-3d electrons 
with spin a in the orbital £ at the site i and Pima- is 
the annihilation operator for As-Ap electrons with spin 
a in the orbital m at the site i. In eq. (2), the transfer 
integrals tf^j^^ g,, if,^- * -.^-f,™ and the atomic energies 
Eg, are determined so as to fit both the energy and 
the weights of orbitals for each band obtained from the 
tight-binding approximation to those from the density 
functional calculation.^'*' The doping x corresponds to 
the number of electrons per unit cell n = 24 + 2a; in the 
present model. 

We show the band structure obtained from the d-p 
tight-binding Hamiltonian eq. (1) together with that ob- 
tained from the density functional calculation in the left 
panel of Fig. 1. It is found that the former reproduces 
the latter very well. We note that the weights of orbitals 
also agree very well with each other (not shown). The 10 
bands near the Fermi level are mainly constructed by the 
Fe 3d orbitals and the 6 bands below the 3d 10 bands are 
mainly constructed by the As 4p orbitals (not shown in 
Fig. 1). The Fermi surface for the d-p tight-binding model 
is shown in the right panel of Fig. 1, where we can see 
nearly circular hole pockets around the F point and ellip- 




Fig. 1. (Color online) (Left panel) The band structure obtained 
from oq. (1) (solid line) and that obtained from the density func- 
tional calculation (open circle). (Right panel) Fermi surface ob- 
tained from the d-p model eq. (1) for x = 
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Fig. 2. (Color online) Several components of the noninteracting 
susceptibility X^'J^j"^^^^ for a; = (a) and those for x = 0.1 (b). 



tical electron pockets around the M point. These results 
are consistent with the first principle calculations. 

Now we consider the effect of the Coulomb inter- 
action on Fe site: the intra-orbital (inter-orbital) direct 
terms U ([/'), the Hund's rule coupling J and the pair- 
transfer J'. Within the RPA,^^' the spin susceptibil- 
ity X*(q) and the charge- orbital susceptibility x'^(q) are 
given in the 50 x 50 matrix representation as follows, 



x<"Hq)(?)-ix^°Hq) 



x^(q) = (1 

with the noninteracting susceptibility 



(2) 
(3) 



(0) Q,/3 



(q) = - 



EE 

,/3 



X u?^,,(k)*<^,^(k + qX ,,(k)<_^,(k + q) , (4) 

where /i, v (=1-16) are band indexes, a, (3 {=A,B) rep- 
resent two Fe sites, £ represents Fe 3d orbitals, M"^(k) is 
the eigenvector which diagonalizes Hq eq. (1), Ck,^ is the 
corresponding eigenenergy of band ji with wave vector k 
and /(e) is the Fermi distribution function. In eqs. (2) 
and (3), the interaction matrix S (C) is given by 



S {C)^ { 



In the weak coupling regime, the superconducting 
gap equation is given by 



{U{U) 
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= /?, 4 = ^2 = h 


= 4) 
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(5) 



where A"^,(k) is the gap function and ^e'lil isiS'i) 
the effective pairing interaction. Within the RPA,^^' 



a, 13 



i\i2 ^3^4^*^) given in the 50 x 50 matrix. 



\s)-\{Cr{^)C-\c\ (6) 



with ^7 = f for the spin-singlet state and 77 : 



for the 



spin-triplet state. The gap equation eq. (5) is solved to 
obtain the gap function A"^'?(k) with the eigenvalue A. 
At T = Tc, the largest eigenvalue A becomes unity. 

In the present paper, we mainly focus on the case 
with X = 0.1, where the superconductivity is observed 
in the compounds.^' For simplicity, we set T — 0.02eV 
and U = U' + 2J , J = J' . We use 32 x 32 k points in 
the numerical calculations for eqs. (2)- (6), and also use 
the first Fourier transformation (FFT) to solve the gap 
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Fig. 4. (Color online) (a) {/'-dependence of the eigenvalue A 
for X = 0.1 and T = 0.02eV. The closed circles and opened 
squares represent A of the s-wave and dxy-'wa.ve, respectively. 
The gap function ^^2^_ 2 2i^) dxy-wavc symmetry at 
U' = 1.77eV,J = OeV and x = 0.1 (b), and that for s-wave 
symmetry at U' = O.SeV, J = 0.3eV and x = 0.1 (c). The solid 
lines represent the nodes of the gap function. 



Fig. 3. (Color online) Several components of the spin suscep- 
tibility Xijj^^i^i^ (3-) s-nd (b), the charge-orbital susceptibil- 
ity Xi^'t!fe^i^ (c) and (d) and the effective pairing interaction 
V"^'tl i,^(^ (e) and (f), for U' = 1.77eV and J = OeV (a) (c) (e), 
and for U' = O.SeV and J = O.SeV (b) (d) (f), respectively. 



equation eq. (5). 

Fig. 2 shows the several components of the nonin- 
teracting susceptibihty given in eq. (4) for a; = 0.1 to- 
gether with those for a; = for comparison. We see the 
peaks centered at the T point and those at the M point 
for X = 0, where the former peaks are due to the nest- 
ing between the hole (electron) pockets and the latter 
peaks are due to the nesting between the hole pockets 
and the electron pockets. With the electron doping, the 
hole pockets around the F point shrink and the smallest 
one disappears for x ~ 0.1, while the electron pockets 
become larger (see Fig. 1 of ref. 23). As the result, the 
nesting effect becomes weak, and then, the peak at the 
M point is suppressed and shifts around the AI point 
for X = 0.1, resulting in the incommensurate SDW as 
mentioned later. 

The several components of the spin susceptibility 
Xii'i^^iii'i) given in eq. (2) are plotted in Figs. 3 (a) 
and (b), where the parameters U' and J are set to the 
condition with A = 1 (see Fig. 4 (a)). The spin suscepti- 
bility is enhanced due to the effect of the Coulomb inter- 
action, especially for the diagonal component of dr^2_.y2. 



For J — OeV, the incommensurate peaks around M point 
are observed as reflecting the structure of the bare sus- 
ceptibility shown in Fig. 2 (b). On the other hand, for 
J = O.SeV, the commensurate peaks centered at the M 
point are observed; which is due to an effect of the Hund's 
coupling J. 

The several components of the charge-orbital suscep- 
tibility Xt^if^i^iS'^) given in eq. (3) are plotted in Figs. 
(3) (c) and (d). In contrast to the case with the spin sus- 
ceptibility, the off-diagonal component of d^2_y2 — dyz 
which corresponds to the orbital susceptibility becomes 
most dominant due to the effect of the inter-orbital 
Coulomb interaction U' . For J = OeV, the orbital sus- 
ceptibility is enhanced and shows peaks around the M 
point together with those at the F point. On the other 
hand, for J — 0.3eV, the enhancement of the orbital sus- 
ceptibility is very small, where U' is smaller than the 
intra-orbital interaction [/ = [/' + 2 J which suppresses 
the orbital susceptibility. 

The several components of the effective pairing in- 
teraction V^^f^ 1^1^ (q) for the spin-singlet state given in 
eq. (6) are plotted in Figs. 3 (e) and (f). Since the 
largest eigenvalue A is always spin-singlet state in the 
present study, we show the effective pairing interaction 
only for the spin-singlet state. Their structures are sim- 
ilar to those of the spin susceptibility in the both cases 
for J = OeV and J = 0.3eV. This is because the contri- 
butions of the spin fluctuations to the effective pairing 
interaction is three times larger than those of the orbital 
fluctuations according to eq. (6). 
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Fig. 5. (Color online) The phase diagram on U'-J plane for x = 
0.1 and T = 0.02eV. The lines represent the instabilities for the 
stripe antiferromagnetic order (solid), the incommensurate SDW 
(dot-dashed), the extended s-wave superconductivity (dotted) 
and the d^y-wave superconductivity (dashed), respectively. The 
inset shows the magnetic phase diagram for a; = 0. 



Substituting ^i^e^ iaeS*^^ S^P equation eq. 

(5), we obtain the gap function A"^'?(k) with the eigen- 
value A. In Fig. 4 (a), the two largest eigenvalues A, which 
are for the extended s-wave and the dxy-wave pairing 
symmetries, are plotted as functions of U' for several 
values of J. We confirmed that eigenvalues for the other 
paring symmetries are much smaller. With increasing U' , 
A monotonically increases and finally becomes unity at a 
critical value U'^ above which the superconducting state 
is realized. For J = OeV, the largest eigenvalue A is for 
the da;y-wave symmetry, where the gap function has line 
nodes on the Fermi surfaces as shown in Fig. (4) (b). On 
the other hand, for J = 0.3eV, the largest eigenvalue A is 
for the extended s-wave symmetry, where the gap func- 
tion changes its sign between the hole pockets and the 
electron pockets without nodes on the Fermi surfaces as 
shown in Fig. (4) (c).",i9,20) p^j. j ^ Q.lSeV, these two 
eigenvalues are almost degenerate. The extended s-wave 
pairing is mediated by the pairing interaction with the 
sharp peak at Q = (7r,7r) (see Fig. 3 (f)), while the dxy- 
wave pairing is mediated by the pairing interaction with 
the incommensurate peaks around Q = (vr, tt) together 
with the peak at Q = (0,0) (see Fig. 3 (e)). 

The phase diagram on U'-J plane for x = 0.1 and 
T = 0.02eV is shown in Fig. 5, where the magnetic in- 
stability is determined by the divergence of the spin sus- 
ceptibility and the superconducting instability is deter- 
mined by A = 1 as mentioned before. For J > 0.25eV, the 
stripe-type antiferromagnetic order with Q = (tt, tt) ap- 
pears, while, for J < 0.25eV, the incommensurate SDW 
(ISDW) with Q ~ (tt, tt) appears (see also Figs. 3 (a) 
and (b)). It is noted that we only observe the stripe- type 
antiferromagnetic order for a; = as shown in the inset 
in Fig. 5. The extended s-wave pairing is realized near 



the stripe-type antiferromagnetic order for J > O.lSeV, 
where the spin fluctuation with Q = (tt, tt) is enhanced 
as shown in Fig. 3 (b). On the other hand, the d^y- 
wave pairing is realized near the ISDW for J < O.lSeV, 
where the spin fluctuation with Q ~ (tt, tt) is enhanced 
as shown in Fig. 3 (a). 

In summary, we investigated the pairing symmetry 
of the two-dimensional 16-band d-p model by using the 
RPA and obtained the phase diagram on the U'-J plane 
for X = 0.1, T = 0.02eV. For J > O.lSeV, the most favor- 
able pairing is extended s-wave symmetry whose order 
parameter changes its sign between the hole pockets and 
the electron pockets, while for J < 0.15eV, it is d^jy-wave 
symmetry. Then, the effect of the Hund's coupling J is 
crucial to realize the extended s-wave pairing in the pure 
system without impurities. According to the recent ex- 
periment of very weak Tc-suppression by Co- impurities, 
we suppose that the d^^y-wave pairing is suppressed by 
pair breaking effect and the extended s-wave paring is 
realized in real materials. 
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